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Abstract. Flavour effects play an important role in the statistical evolution of particle number 
densities in several particle physics phenomena. We present a fully flavour-covariant formalism 
for transport phenomena, in order to consistently capture all flavour effects in the system. We 
explicitly study the scenario of Resonant Leptogenesis (RL), and show that flavour covariance 
requires one to consider generically off-diagonal number densities, rank-4 rate tensors in flavour 
space, and non-trivial generalization of the discrete symmetries C, P and T. The flavour- 
covariant transport equations, obtained in our semi-classical framework, describe the effects of 
three relevant physical phenomena: coherent heavy-neutrino oscillations, quantum decoherence 
in the charged-lepton sector, and resonant CP violation due to heavy-neutrino mixing. We show 
quantitatively that the final asymmetry predicted in RL models may vary by as much as an order 
of magnitude between partially flavour off-diagonal treatments. A full field-theoretic treatment 
in the weakly-resonant regime, based on the Kadanoff-Baym (KB) equations, confirms that 
heavy-neutrino oscillations and mixing are two distinct phenomena, and reproduces the results 
obtained in our semi-classical framework. Finally, we show that the quasi-particle ansaetze, 
often employed in KB approaches to RL, discard the phenomenon of mixing, capturing only 
oscillations and leading to an underestimate of the final asymmetry by a factor of order 2. 


1. Introduction 

Leptogenesis [I] is an elegant unifying framework for dynamically generating both the measured 
matter-antimatter asymmetry in our Universe and the observed smallness of the light neutrino 
masses |2|. In scenarios of Resonant Leptogenesis (RL) |3j|4|, this mechanism may be testable 
in foreseeable laboratory experiments. RL relies on the fact that the e-type CP-asymmetry 
becomes dominant [ 5 ] and gets resonantly enhanced, when at least two of the heavy neutrinos 
have a small mass difference comparable to their decay widths |3|. This resonant enhancement 
allows a successful low-scale leptogenesis |4j|61, whilst retaining perfect agreement with the light- 
neutrino oscillation data. The level of testability is further extended in the scenario of Resonant 
CGenesis (RL^), where the final lepton asymmetry is dominantly generated and stored in a single 
lepton flavour t [7j[8j. In such models, the heavy neutrinos could be as light as the electroweak 
scale [6j, whilst still having sizable couplings to other charged-lepton flavours i' / l. Thus, RL^ 
scenarios may be directly testable at the energy frontier in the run-II phase of the LHC [9j, as 



well as in various low-energy experiments searching for lepton flavour/number violation 110 at 
the intensity frontier. 

Flavour effects in both heavy-neutrino and charged-lepton sectors, as well as the interplay 
between them, can play an important role in determining the final lepton asymmetry in low- 
scale leptogenesis models (for a review, see e.g. 111]). These intrinsically-quantum effects can, 
in principle, be accounted for by extending the classical flavour-diagonal Boltzmann equations 
for the number densities of individual flavour species to a semi-classical evolution equation for 
a matrix of number densities , analogous to the formalism presented in 12 for light neutrinos. 


This so-called ‘density matrix’ formalism has been adopted to describe flavour effects in various 
leptogenesis scenarios 113 15]. It was recently shown [l6|, in a semi-classical approach, that 
a consistent treatment of all pertinent flavour effects, including flavour mixing, oscillations 
and off-diagonal (de)coherences, necessitates a fully flavour-covariant formalism, in order to 
provide a complete and unified description of RL; for a summary, see |17|. In this flavour- 
covariant formalism, the resonant mixing of different heavy-neutrino flavours and coherent 
oscillations between them are found to be two distinct physical phenomena, in analogy with 
the experimentally-distinguishable phenomena of mixing and oscillations in the neutral K-, D-, 
B- and £> s -meson systems [2], 

One can go beyond the semi-classical ‘density-matrix’ approach to leptogenesis by means 
of a quantum field-theoretic analogue of the Boltzmann equations, known as the Kadanoff- 
Baym (KB) equations |l8] (for a review, see e.g. (19)1). Such ‘first-principles’ approaches to 
leptogenesis 1201 are, in principle, capable of accounting consistently for all flavour effects, 
in addition to off-shell and finite-width effects, including thermal corrections. However, it is 
often necessary to use truncated gradient expansions and quasi-particle ansaetze to relate the 
propagators appearing in the KB equations to particle number densities. Recently, using the 
novel perturbative formulation of thermal field theory developed in 121 , it was shown 22 that 


quantum transport equations for leptogenesis can be obtained from the KB formalism without 
the need for gradient expansion or quasi-particle ansaetze, thereby capturing fully the pertinent 
flavour effects. Specifically, the source term for the lepton asymmetry obtained, at leading order, 


in this KB approach 22 was found to be exactly the same as that obtained in the semi-classical 
flavour-covariant approach of 16 , confirming that flavour mixing and oscillations are indeed 
two physically-distinct phenomena. The proper treatment of these flavour effects may have a 
significant effect upon the final lepton asymmetry, as compared to partially flavour-dependent 
or flavour-diagonal limits, thereby altering the viable parameter space for models of RL and 
impacting upon the prospects of testing the leptogenesis mechanism. 

The plan of these proceedings is as follows. In Section [2j we review the main features of 
our fully flavour-covariant formalism in the context of leptogenesis. In Section [3j we present the 
Markovian flavour-covariant transport equations for lepton and heavy-neutrino number densities 
with arbitrary flavour content. We also discuss a numerical example to illustrate the full impact 
of the flavour off-diagonal effects within the context of an RL r model. In Section]!] we derive the 
quantum transport equations relevant to the source term for the lepton asymmetry, following a 
well-defined perturbative loopwise truncation scheme, making comparison with the semi-classical 
approach discussed in Section]!] Our conclusions are given in Section]!] 


2. Flavour-covariant formalism 

We consider the lepton-doublet field operators Li (with l = 1,2,...,A f£) and right-handed 
Majorana neutrino field operators 1 Vr ;Q (with a = 1,2,... ,AI’m), with arbitrary flavour content, 
transforming as follows in the fundamental representation of U(Nl) <S> U(AfN): 
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where V) m 6 U(Nl) and U a ^ £ U(Afi\r). The relevant neutrino Lagrangian is given by 

~C N = V^ JV R,a + ^,aN% + H.C . 1 

where $ = *(72$* is the isospin conjugate of the Higgs doublet $ and the superscript C denotes 
charge conjugation. The Lagrangian © transforms covariantly under U(J\Tl)®U(J\Tn), provided 
the Yukawa couplings and Majorana mass matrix transform as 

K h\ a = vru a phj, [M n }<*P -> [M' n ]^ = u% u fi s [M N y s . (3) 

In this flavour-covariant formalism, the plane-wave decompositions of the field operators are 
written in a manifestly flavour-covariant way, e.g. 
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( 4 ) 


where f p = f s is the helicity index and [E\( p)] ; m = p 2 (5 z m +[M^Mi]j m , with Ml being the 
charged-lepton mass matrix. Thus, flavour covariance requires the Dirac four-spinors [u(p, s)] 
and [u(p, s)\j k to transform as rank-2 tensors in flavour space. The creation and annihilation 
operators b k = b\, bk, df’ k = dk and d\ satisfy the equal-time anti-commutation relations 

{bi(p,s,t), b m (p',s',t)} = {d^ m (p,s,t), d](p',s',t)} = (27r) 3 <5 (3) (p - p') S 3s > . (5) 

For the heavy Majorana neutrino creation and annihilation operators a“(k, r, t) and a Q (k, r, t), 
it is necessary to introduce the flavour-covariant Majorana constraint 

rf t,Q (k, — r,t) = G a P bp(k,r,i) = G Q/3 a /3 (k,r,t) , (6) 

where G a0 = [U*U^] aP are the elements of a unitary matrix G, which transforms as a 
contravariant rank-2 tensor under f/(A/jv)- Notice that G = 1 in the mass eigenbasis. 

Similar flavour rotations are forced by flavour-covariance under the discrete symmetry 
transformations C, P and T. This necessarily leads to generalized C and T transformations: 


bi(p,s,t) 5 = G lm b m (p,s,t) c = -id^ l (p,s,t) , (7a) 

bi(p,s,t) p = -sbi(-p, -s,t) , (7b) 

bi(p,s,t) T = Gim b m (p, s, i) T = bi(-p,s,-t) , (7c) 


where G lm = [V*V^] lm is the lepton analogue of the heavy-neutrino tensor G. 

We may now define the matrix number densities of the leptons and heavy neutrinos, which 
describe completely the flavour content of the system: 

Ks 2 (P,t)\r = V 3 _1 ( &m (p ,S 2 ,t)bl{p,S 1 ,t))t , 

[nf 1 S 2 (p,t)]r = V^ 1 {d}(p,s 1 ,t)d^ m {p,s 2 ,t))t , 

[ n nr 2 (k; t)\cf = V 3 -1 (a P (k, r 2 , t)a a (k, n, t)) t . 


(8a) 

(8b) 

(8c) 





Here, V 3 is the infinite coordinate three-volume and t = t — t t is the macroscopic time, equal 
to the time interval between specification of initial conditions (U) and subsequent observation 
of the system (t). The total number densities n x (t) are obtained by tracing over helicity and 
isospin indices and integrating over the three-momenta. For the Majorana neutrinos, n N and 
n N are not independent quantities and are related by the generalized Majorana constraint <©• 
Using the generalized discrete transformations (7a)-(7c), we can define the flavour-covariant 
CP-“even” and-“odd” quantities 
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Sn N = n N — n N , 8n L = n L — n L , (9) 

which will be used in the following section to write down the flavour-covariant rate equations. 


3. Flavour-covariant semi-classical rate equations 

We first derive a master equation governing the time evolution of the matrix number densities 


n x (p,t), as given in ( 8 a)—( 8 c). By using the Liouville-von Neumann and Heisenberg equations 


of motion and subsequently performing a Wigner-Weisskopf approximation in the Markovian 
limit 1161 , we find 


d , 1 

—n x (k, t) ~ n x (k,t)]) t - -J dt' ([H int (t'), n x {k,t)]]) t , (10) 

where H x and Pint are respectively the free and interaction parts of the Hamiltonian and 

n x (k,t) = (n x (k,t-ti)) t = Tr{p(t;ti)h x (k,t;ti)} (11) 

is the ensemble expectation value of the quantum-mechanical number-density operator 
fi x (k, t: ti), in which p(t: ti) is the interaction-picture density operator. The first term on the 
RHS of ( |10[ ), involving the free Hamiltonian, is responsible for flavour oscillations, whereas the 
second term contains the collision terms of the generalized Boltzmann equations. Explicitly, for 


the system of charged leptons and heavy neutrinos, we find 16 
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( 12 a) 


^ [< r2 (k,t)]/ = -i[E N { k), n x r2 (k,t)]\ [C x r2 (k,t)]J + G a \ [< ri (k,t)]/G^ . (12b) 


The collision terms [C 8 s {p, t)]™ and [C^ r (k, t)]a involve the product of two new rank-4 
tensors in flavour space, namely, the statistical number density tensor and the absorptive rate 
tensor, whose appearance is necessary for the flavour covariance of the formalism 16 . The 


emergence of these rank-4 tensors may be understood in terms of the unitarity cuts of the 
partial self-energies, as was shown by an explicit calculation of the relevant transition amplitudes 
using a generalized optical theorem in 16 . The off-diagonal components of the rate tensor are 


responsible for the evolution of flavour-coherences in the system. 

In the limit when two (or more) of the heavy Majorana neutrinos become degenerate, the 
e-type CP- violation can be resonantly enhanced, even up to order one [3|, due to the interference 
between the tree-level and self-energy-corrected decays. In this regime, finite-order perturbation 
theory breaks down and one must resum the self-energy corrections in order to account for the 
heavy-neutrino mixing effects. In the semi-classical approach, we perform such resummation in 
an effective way by replacing the tree-level neutrino Yukawa couplings h l by their resummed 
counterparts h ^ in the transport equations. In the next section, this approach will be justified 













for the weakly-resonant regime of RL, where Ihv a(3 < \ m N a — m/v a m N a p , by using a 
‘first-principles’ field-theoretic approach. The explicit algebraic form of the resummed neutrino 
Yukawa couplings in the heavy-neutrino mass eigenbasis can be found in |4 ; the corresponding 

form in a general flavour basis may be obtained by the appropriate flavour transformation, 

P — 


i.e. h ; “ = V l m U a phrrf, where h m p = h m p in the mass eigenbasis 

In addition, we make the following reasonable approximations to simplify the flavour-covariant 
rate equations for RL: we (i) assume kinetic equilibrium, which is ensured by the presence of fast 
elastic-scattering processes; (ii) work in the classical-statistical regime; (iii) neglect thermal and 
chemical-potential effects |6j; and (iv) neglect the mass splitting between different heavy-neutrino 
flavours inside thermal integrals, using an average mass mjv and energy Efy(k) = (|k| 2 +m^) 1//2 , 
as is appropriate since the average momentum scale |k| ~ T \niN a ~ mjv a |. In order to 
guarantee the correct equilibrium behaviour, we must also include the effect of the A L = 0 and 
A L = 2 scattering processes, with proper real intermediate state (RIS) subtraction [4 16 23 . As 
detailed in |16] , it is necessary to account for thermal corrections in the RIS contributions, when 
considering off-diagonal heavy-neutrino flavour correlations. Finally, it is important to include 
the effect of the charged-lepton Yukawa couplings, which are responsible for the decoherence of 
the charged leptons towards their would-be mass eigenbasis. 

Taking into account the expansion of the Universe, we then obtain the following manifestly 
flavour-covariant rate equations for the normalized CP-“even” number density matrix rj N and 
CP- “odd” number density matrices Srj N and Sr] L (where r] X = n x /n 7 , with n 7 being the photon 
number density) fit 


16 



Here, is the Hubble parameter at z = m^/T = 1 and we have defined the thermally-averaged 
heavy-neutrino energy matrix 


Sn = 



(14) 


where = 2 counts the helicity degrees of freedom. In addition, and are respectively 
the CP-“even” and -“odd” thermally-averaged rate tensors, describing heavy-neutrino decays 
and inverse decays and written in terms of the resununed Yukawa couplings as 


hUJ = f (hah 15 + ruirf) . IS-&IJ = f (hp - [h E ]„[h 5 ]' 5 ) , (15) 
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where c denotes the generalized CP-conjugate and we have used the shorthand notation 
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with the phase-space measure for the species X given by 
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In (13c), and respectively describe the washout due to A L = 0 and A L = 2 

resonant scattering, and 7d ec and <57dec k govern the charged-lepton decoherence 16 . Lastly, 


for a Hermitian matrix A, we have defined the flavour-covariant generalized real and imaginary 
parts 


[Re(A)]/ = ^(V + G^Vg-O), 
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-G aX A*G^y (18) 


which reduce to the usual real and imaginary parts in the heavy-neutrino mass eigenbasis. 


The flavour-covariant rate equations (13a)—(13c) provide a complete and unified description 


of the generation of the lepton asymmetry in RL, consistently describing the following physically- 
distinct effects in a single framework: 

(i) Lepton asymmetry due to the CP-violating resonant mixing between heavy neutrinos, as 
described by the CP- “odd” rates S'y ^, appearing in the first two terms on the RHS of (13c). 


This provides a flavour-covariant generalization of the mixing effects discussed earlier in |4j. 
(ii) Lepton asymmetry via coherent heavy-neutrino oscillations. This is an 0(/i 4 ) effect on the 

, thus differing from the (D(h e ) mechanism studied in 


total lepton asymmetry 


16 


15 


which 


typically takes place at temperatures much higher than the sterile neutrino masses. 

(iii) Decoherence effects due to charged-lepton Yukawa couplings. Our description of these effects 
generalizes the analysis in 1131 to an arbitrary flavour basis. 


In order to illustrate the importance of the flavour effects captured only by the flavour- 
covariant rate equations ( 13a)~( [T3c| ), we consider an RL r model, comprising an approximately 
•SO^-symmetric heavy-neutrino sector at the grand unification scale px ~ 2 x 10 16 GeV, with 
masses Mx(px) = m/vl 3 + AM/v |6|. The soft 5'0(3)-breaking mass term is taken to be of 
the form A= diag (AMi, AM 2 / 2 , — AM 2 / 2 ). At the electroweak scale, an additional mass 
splitting arises from the RG running, such that lL/jv(mjv) = mArl _3 + A+ AIn order 
to accommodate the smallness of the light neutrino masses in a technically natural manner, 
we also require the heavy-neutrino Yukawa sector to possess an approximate leptonic 17(1)/ 
symmetry. This results in the following structure for the heavy-neutrino Yukawa couplings: 
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(19) 


where a,b,c are arbitrary complex parameters. If the theory were to have an exact 17(1)/ 
symmetry, i.e. if the 17(l)/-breaking parameters e e ,n,T were set to zero, the light neutrinos would 
remain massless to all orders in perturbation theory. For electroweak-scale heavy neutrinos, we 
require |a|, |6| < 10 -2 , in order to be consistent with current light-neutrino mass bounds, and 


|c| < 10 -5 and \e e ^ )T \ < 10 , in order to protect the r lepton asymmetry from wash-out effects. 

Using the Yukawa coupling given by (19), we solve the rate equations (13a)-(13c) numerically 
to obtain the total lepton asymmetry 5rp = Tr(Sri L ) in our flavour-covariant formalism. This 
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-3 x 10~ 5 
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(8.04 -3.79 i) x 10" 3 

e e 

5.73 * x 10~ 8 


4.30* x 10~ 7 


6.39* x 10~ 7 
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Figure 1: (a) The total lepton asymmetry for a minimal RL r model with the parameters in (b). 
The thick solid lines show the total asymmetry obtained in our fully flavour-covariant formalism 
for different intial conditions; the dashed lines show that obtained in various flavour-diagonal 
limits. 


as given in Figure lb 
Figure 


la 


is shown in Figure la for a typical set of benchmark values for the Yukawa coupling parameters, 

which is consistent with all current experimental constraints [l6]. In 
the horizontal dotted line shows the value of 5r] L required to explain the observed 
baryon asymmetry, whereas the vertical line shows the critical temperature z c = mjsr/T c , beyond 
which the electroweak sphaleron processes become ineffective in converting lepton asymmetry 
to baryon asymmetry. The thick solid lines show the evolution of 5i] L for three different 
initial conditions, to which the final lepton asymmetry 5r] L (z 1) is insensitive as a general 
consequence of the RL mechanism in the strong-washout regime jb|. For comparison, Figure [Ia| 
also shows various flavour-diagonal limits, i.e. when either the heavy-neutrino (dashed line) or 
the lepton (dash-dotted line) number density or both (dotted line) are diagonal in flavour space. 
Also shown (thin solid line) is the approximate analytic solution obtained in |16| for the case 
of a diagonal heavy-neutrino number density. The enhancement of the lepton asymmetry in 
the fully flavour-covariant formalism (solid line), as compared to assuming a flavour-diagonal 
heavy-neutrino number density (dashed line), is mainly due to coherent oscillations between 
the heavy-neutrino flavours, leading to a factor of 2 increase. Finally, we observe that the 
predicted lepton asymmetry differs by approximately an order of magnitude between the two 
partially flavour off-diagonal treatments (dashed and dash-dotted lines). This provides a striking 
illustration of the importance of capturing all pertinent flavour effects and their interplay by 
means of a fully flavour-covariant formulation for transport phenomena. 


4. Flavour-covariant quantum transport equations 

The semi-classical approach detailed in the preceding section has the advantage that it is 
constructed with physical observables, i.e. particle number densities, in mind. However, it has 
the disadvantage that quantum effects, such as finite particle widths, must be incorporated in 
an effective manner. For instance, in the latter example, we must subtract the RIS contributions 


from the collision terms 23 , which would otherwise lead to double counting of decay and inverse- 





















decay processes. It is therefore desirable to seek a more first-principles description of transport 
phenomena, in which quantum effects are incorporated consistently from the out-set. 


Such a description is provided by the Kadanoff-Baym (KB) approach [18] (see also 19,24] 


constructed within the Schwinger-Keldysh closed-time path (CTP) formalism of thermal field 
theory [251. Therein, one arrives at systems of KB equations by partially inverting the Schwinger- 
Dyson equation of the 2PI CJT effective action 26 . Unfortunately, the KB equations describe 


the spacetime evolution of propagators and, as a result, it is necessary to use approximation 
schemes in order to obtain the quantum transport equations of particle number densities. 

In order to avoid the technical complications of spinor fields, we will consider a scalar model 
of RL, see 22 , comprising: two real scalar fields N a (cc = 1,2), modelling heavy-neutrinos of 


two flavours; one complex scalar field L, modelling charged-leptons of a single flavour; and a 
real scalar held <h, modelling the Standard Model Higgs. The lepton number can be associated 
with the global 1/(1) symmetry of the complex scalar held L. 

The CTP formalism may be formulated in two ways, working either in the Heisenberg or 
interaction picture. In the former (see e.g. 27 ), the density operator p(tf,ti) does not evolve 
in time, remaining hxed at the boundary time t*. As a result, the free propagators encode the 
initial conditions of the statistical ensemble, e.g. 


[*A>’°(a:,y,0)] ( 
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= (N a (x)Nf*(y)) 0 = -TV p(0)N a (x)N^(y) . 
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( 20 ) 


In this case, it is well-known that there does not exist a well-defined perturbative expansion. This 
may be understood by considering the Taylor expansion of the exponential decay to equilibrium: 
e _ri = 1 —rt+^(rt) 2 -|----. Any truncation of this expansion at a finite order in the decay rate 
r leads to secular behaviour when t > 1/T 1191. This problem manifests in the Feynman-Dyson 
series of the Heisenberg interpretation of the CTP formalism as pinch singularities 28 , which 
result from ill-defined products of delta functions with identical arguments. As a consequence, 
it is necessary to work with dressed propagators, and therefore, ansaetze are required in order to 
extract particle number densities. The most common is the quasi-particle approximation known 
as the KB ansatz: 


N A < (M)]/ = 2vr5(A: 2 - ?n^)[n^ ress (k,t)] Q /3 , k° > 0 , (21) 

where [ra^ ess (k, t)]a? approximates the matrix number density of spectrally-dressed particles. 

On the other hand and in strong contrast to the earlier literature, it was shown recently 121 
that a perturbative framework of non-equilibrium thermal field theory is in fact viable, if we 
work instead in the interaction picture. Since the interaction-picture density operator p(tf,ti ) 
evolves in time, being evaluated at a macroscopic time t = tf — ti after the specification of the 
initial conditions, the free positive-frequency Wightman propagator becomes 

[iA>'°(x,y,tf,ti)\ a 0 = ( N a(x-,t i )N l3 (y;ti)) t = TV p(t f ; ii)N a (x\U)N^(y; U) . (22) 

In the heavy-neutrino mass eigenbasis and assuming spatial homogeneity, the free Wightman 
propagators then have the following explicit form: 

[iA^’°(fc, k', tf, ti)\a/3 = 2Tr\2ko\ 1 ^ 2 5(k 2 — mjv, a ) 27 t|2A;q| 1 / 2 <5(&: /2 — rriN,p) e^ ko ~ k °^ f 
x (d{±k 0 )e(±k' Q )5 af} + [d(fc o )0(A:( ) ) + 0(-A:o)0(-feo)][n iV (k,t)] Q/3 )(27r) 3 ( 5( 3 )(k-k / ) , (23) 

depending on the time-dependent matrix number density n N (k, t) of spectrally-free particles. 
These number densities appear as unknown functions in the Feynman-Dyson series, with their 








functional form being fixed only after the governing transport equations have been solved. Thus, 
the exponential decay to equilibrium is present implicitly in the free propagators of the theory, 
thereby avoiding the problem of secularity or pinch singularities, see [21 . 

Making a Markovian approximation and additionally setting tun, a — m N = (^ jv,i + wjv , 2 )/ 2 , 
valid in the weakly-resonant regime, the free Wightman propagators in ( |23| ) reduce to 

[iA*'°(k,t)]J = 27r5(k 2 - m 2 N )(d(±k 0 )dJ + [n N (k,t)] a ^ , (24) 


written here in a single-momentum representation and in a general flavour basis. The Markovian 


and homogeneous form of the free propagator in (24) should be compared (for k o > 0) with the 


KB ansatz of the dressed propagator in (21), wherein we note that their spectral structure is 


identical in spite of the fact that the latter should be fully dressed spectrally. 

In coordinate space, the KB equations take the following generic form (see e.g. (24)): 
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a^ = — — (n> * a< — n< * a> + 2 n> * a 
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(25) 


where = d x nd x is the d’Alembertian operator and • indicates matrix multiplication in flavour 
space. The * denotes the convolution 


A * B — 


'zGflt 


A(x, z, tfji ) • B(z, y, tf,ti ) , 


(26) 


which is performed over the hypervolume fit = [U,tf \ xt - [—|, |] x R , bounded temporally 


from below and above by the boundary and observation times, respectively 21 


The KB equation (|25|) may be recast in a double momentum-space representation as follows: 

(27) 


• + xi-p ^ — — — ^n> * — n< ~k + 2 'k A-p 

Here, * denotes the weighted convolution integral 

A* B = [ (2tt ) 4 6^\q-q') A(p,q,tf;U) ■ B(q',p',tf,ti) , 
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Stiqo-q'o) - I Sin ^o)*/2] 


(30) 

7T qo - q 0 

By considering the Noether charge, see [2l], the number density n(f, X) may be related 
unambiguously to the negative-frequency Wightman propagator via 


/•(Aj _ __ r(Z.) r 

n(t, X) = / (p 0 +p' 0 )iA < (p,p',t f -ti) , / = / e ~ l (p-p')' x e(p 0 + j/ 0 ) , (31) 

J p, p' J p, p' J p, p' 

with f p = f and f p p/ = f p f p/ . Following [21], we may then translate (27) into the final rate 


jp,p r jp jp 

equation for the number density 
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1 Z-W/ N 

-J ({n>, a<}* - {n<, a>}* + 2[n<, a p ]*J , 


df 


/p,p' 


(32) 










where we use a compact notation 


[. A , £>]* = A-k B - B ★ A , {A, £?}* = A * T? + -B * A 


(33) 


The first two terms on the LHS of (32) are the drift terms and the latter two account for mean- 


field effects, including oscillations; the terms on the RHS describe collisions. It should be stressed 


that (32), obtained without employing a gradient expansion or quasi-particle ansatz, is valid to 


any order in perturbation theory and accounts fully for spatial inhomogeneity, non-Markovian 
dynamics (memory effects) and flavour effects. 

As identified in 1211, the general rate equation in (32) may be truncated in a perturbative 


loopwise sense in two ways: (i) spectrally: by truncating the external leg, we determine what is 
being counted, e.g. inserting free propagators, we count spectrally-free particles; (ii) statistically: 
by truncating the self-energies, we determine the set of processes that drive the statistical 
evolution, e.g. inserting one-loop self-energies, we include decay and inverse-decay processes. 

4-1. Heavy-neutrino rate equations 

Assuming spatial homogeneity and absorbing the principal part self-energy ITp into the thermal 
mass = |rri 7 v| 2 —lip, we find the rate equation of the dressed heavy-neutrino number density 

-P0 


dn dress 

d t 


f 

Jk, 


k' 


~i[M\ 


2 n , ^ ] - 


1 


in* iA*}, - {in* iA^} 


(34) 


Herein, we have also neglected the commutator involving Ap on the RHS of (|32|), since, in the 


weakly-resonant regime, it contains higher-order effects that are not relevant to this analysis. 

Neglecting 0(h e ) terms proportional to the lepton asymmetry, it is sufficient to approximate 
the charged-lepton and Higgs propagators, appearing in the heavy-neutrino self-energies, by 
their quasi-particle (narrow-width) equilibrium forms with vanishing chemical potential, i.e. 

iA|’ eq (g) = 2n6(q 2 - Ml) [9(±q 0 ) + ra^(q))] , (35) 

iA^ eq (p) = 2n5(p 2 - Ml) [9(±p 0 ) + 0(p o )n^ q (p) + 9(-p 0 )n^ q (p)] • ( 36 ) 

Here, n^(p) = ()eV p2 + M x/' r _ x)^ 1 is the Bose-Einstein distribution and Mx is the thermal 
mass of species X. We are then left with the non-Markovian heavy-neutrino self-energies 

\iU.^{k,k\i r Ma = 2R e(rfh)J 

x [ (2 TT) A 6 t (k -p-q) (2vr )%(k" - p - q) A < eq (p) A5’ eq (g) . 


(37) 




We now perform a Wigner-Weisskopf approximation by replacing f \ by Hoo in all spacetime 
integrals. In the double-momentum representation, this amounts to using the limit 


lim S t (k 0 -po- go) = S(k 0 - p 0 - go) ■ 


t—> OO 


(38) 


At the same time, we replace e fc o)*/ A*(/c, k', tf, ti) —> A<(k, k r , t), absorbing the free- 

phase evolution, which cancels that in the measure of © in the energy-conserving limit. We 
then arrive at the Markovian rate equation for the dressed heavy-neutrino number density 
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(39) 













where the Markovian self-energies 


= 2Re(/ i t h)jB^(k) (40) 

may be written in terms of the thermal loop functions 

B e ^(k) = [ (2ir) 4 8^(p-k + q)Al’ e \q)A^(p) , B^(-k) = B e *(k) € R . (41) 

Subsequently, in the classical-statistical limit, these thermal loop functions can be written as 
P> q (fc 0 > 0, k) = - I dlhj, I dn L (2vr) 4 S^(k - p* - p L ) , (42) 

B^(ko > 0, k) = - J dll* J dn L (2vr) 4 5 (4 \k - p* - p L ) n* (£7*) <(P L ) . (43) 


Since we are interested in the asymmetry at 0(h 4 ), we may truncate (39) spectrally at 


zeroth loop order, replacing the external heavy-neutrino propagators by the free homogeneous 
propagator in (24). The k! integral in (39) can then be performed trivially and we obtain the 


following rate equation for the spectrally-free number density [ n N ] a ®'■ 

- ( ({>n*(k), »A- {in"(*), iA*"(k,()}/)} . (44) 


,N] P 
Jq 

d t 


After substituting for the explicit form of the free heavy-neutrino propagator given by (|24|), we 

i], giving 


assume kinetic equilibrium, as described in 16 , giving 
d [n N ]cf 


d t 




£n, n 
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+ [Re( 7 iT j )]^ - 2^v {"*> Re (7ii (0) )} c • (45) 

eq 


Here, we have defined the CP-“even” rate = /nl$> 2 k a h^, where h a are the tree- 

level Yukawa couplings. In addition, P\- is the thermally-averaged effective energy matrix, given 


by (14) evaluated with thermal masses and gN = 1- 


Lastly, we separate the CP-“even” and - “odd” parts of (|45|), giving 

{n. N , Re(7S° ) )} Q /3 , (46a) 
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= - 2 i 


£n, a 
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1 0 




(46b) 


These final heavy-neutrino rate equations agree, up to 0(h 4 ) as considered here, with those 


obtained by the semi-classical approach [cf. (13a)—(13b)]. 


4-2. Source term for the asymmetry 

In this section, we describe the explicit form of the dressed negative-frequency heavy-neutrino 


Wightman propagator, as derived in 22 using the perturbative approach of 21 , making 


comparison with the standard quasi-particle or KB ansatz. In so doing, we will illustrate 



























Figure 2: Diagrammatic representation of the dressed negative-frequency heavy-neutrino matrix 
Wightman propagator zA<. Double lines indicate fully-dressed propagators and single lines, 
propagators dressed with only dispersive corrections. Unshaded blobs are the relevant self¬ 
energies and shaded blobs are the amputated self-energy corrections to the vertices. 


that both heavy-neutrino mixing and oscillations provide distinct contributions to the 0 (h 4 ) 
lepton asymmetry in the weakly-resonant regime and that the contribution of flavour mixing is 
discarded when the standard quasi-particle approximations are employed. 

Working again in the Markovian regime and assuming that the charged-lepton and Higgs 
propagators have the equilibrium forms in (35) and (36), the dressed negative-frequency heavy- 
neutrino Wightman propagator is determined by the following Schwinger-Dyson equation 21 : 

iA*(k,tf,t) = iA*’°(k,k',t) + iA*’°(k) • in<(fc)( 27 r) 4 h ( 4 ) (A: - k') ■ i A^(fc') 

+ iA^’°(k)-iU R (k)-iA^(k,k',t) + iA<’°(fc, k', t) ■ iIl A (k') ■ iA%(k') . (47) 

In addition, the Schwinger-Dyson equation for the advanced propagator has the closed form 

iA A (k) = iA^°(k) + iA^°(k) • iU A (k) • zA^(fc) . (48) 


As shown in [22], (47) and (48) can be solved iteratively to give an explicit form for the dressed 


negative-frequency Wightman propagator, shown diagrammatically in Figure [2] 


+ 


OO n 

E [0'Ato-in^(fc)) m ] a 7 [zA^°(fe,fc',t)] 7 5 E [«(fc')-<-°(fc')) n ]; • (49) 


m = 0 


n = 0 


The first term on the RHS of (49) contributes to the washout from A L = 0 and A L = 2 
scatterings. The second term instead contributes to the source term for the asymmetry. Rotating 
first to the mass eigenbasis, it was shown in 
of the resummed Yukawa couplings h Q [4,16 


by virtue of the equivalence 


L n = 0 


E ?:A h(* 0 • *n R (fe 


22 ] that this contribution may be written in terms 

(50) 




In addition, the same contribution can be recast as 

[iA= [Ag(fc)]„ 7 ([A^ 0 (fc)]^[iA?’»(t,fc',*)] 1 j[Ai'’ 0 (fc')] j y)[Aj'(t')] w . (51) 


On the other hand, the flavour-covariant KB ansatz for the heavy-neutrino propagator is 

k',t )\ af 3 = 2n5(p 2 - Mjf a )2Tv5(k ' 2 - M^ JS )[n^ B (k,t)] a 0 (2n) 3 5 { 3 ) (k-'k') , (52) 
where we have restricted ko > 0. It is then clear that the KB ansatz satisfies 

(Ar — MjE) [*Ak Bj < {k, k', t )\ a0 = 0, [iA^ B< (k,k',t)] a 0 {k ' 2 — Al‘j i ^ 0 ) = 0, (53) 





















Figure 3: Diagrammatic representation of the factorization of absorptive transitions into the 
resummed Yukawa couplings in the source term for the charged-lepton asymmetry. 


whereas the dressed heavy-neutrino Wightman propagator in (49) does not, due to the mixing 


that gives rise to the resummed Yukawas. We therefore conclude that the KB ansatz discards 
the phenomenon of flavour mixing, accounting only for the separate phenomenon of oscillations. 


Using the equivalence in (50), the effect of mixing, arising from absorptive transitions in the 


dressed propagator (49), can be absorbed into the resummed Yukawa couplings. Once this is 


done, it is then appropriate to replace the non-homogeneous free heavy-neutrino propagator 
A<’°(fc, k', t), appearing on the RHS of (49), by its homogeneous counterpart, as given by (24). 
This procedure is illustrated diagrammatically in Figure [3j We then find that the contribution 
from the charged-lepton self-energy to the source term for the asymmetry is given by 


d bn 1 
d t 


D - 


^(A:o)[h^h“([*A^ 0 (A : ,t)]/^(fc)-[iA^ 0 (fc,t)]/5^(/ ! )) - C.c] . (54) 


We again assume kinetic equilibrium and separate out the CP-“even” and “-odd” parts of the 
heavy-neutrino number density, i.e. n N and Sn N , giving the final form 


d5n L 
d t 


(I="K - sA [<•&]„“ + [*]/ + I, (55) 

V '^eq / z '^eq 


where W[Sn L ] denotes washout terms not studied explicitly here. The source term in (55) agrees 


with (13c), as derived in the semi-classical approach in [16|, up to 0(/i 4 ). Evidently, both flavour 
mixing and oscillations can be identified in (55) as providing distinct contributions to the final 
asymmetry. 

The self-consistent loopwise perturbative truncation scheme described in this section is 
summarized diagrammatically in Figure [4j Before concluding, however, it is important to stress 
the following. In the weakly-resonant regime, while it is appropriate to truncate the heavy- 
neutrino rate equations spectrally at zeroth order, inserting free heavy-neutrino propagators in 
the external legs, it is not appropriate to insert the same free heavy-neutrino propagators in 
the charged-lepton self-energies of the source term for the asymmetry. The latter would instead 
correspond to a zeroth-order statistical truncation, which, were it performed, would have the 
same impact as making a KB ansatz, discarding the physical phenomenon of flavour mixing. 


5. Conclusions 

We have presented a novel approach to the study of flavour effects in RL by means of a fully 
flavour-covariant formalism for transport phenomena 16 . Our manifestly flavour-covariant 


rate equations for heavy-neutrino and lepton number densities provide a complete and unified 
description of RL, capturing three distinct physical effects: (i) resonant mixing between the 
heavy-neutrino flavours, (ii) coherent oscillations between different heavy-neutrino flavours and 
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Figure 4: Diagrammatic representation of the truncation procedure, proceeding spectrally for the 
heavy-neutrino rate equation (top) and statistically for the charged-lepton asymmetry (bottom). 


(iii) quantum decoherence effects in the charged-lepton sector. The full impact of these flavour 
off-diagonal effects has been illustrated in an RL r model, where the total lepton asymmetry 
varies as much as an order of magnitude between the partially flavour off-diagonal treatments. 

We have also presented an embedding of our flavour-covariant formalism 116] within a 
perturbative formulation of non-equilibrium thermal field theory 121 , enabling us to extract 
physically-meaningful particle number densities at any order in perturbation theory, without 
the need for quasi-particle ansaetze. In this novel quantum field-theoretic approach to the KB 
formalism 122], we have justified, at leading order and in the weakly-resonant regime, the semi- 
classical Boltzmann approach adopted in |16], capturing all flavour effects pertinent to RL. In 
particular, we have confirmed that the mixing and oscillations between different heavy-neutrino 
flavours are two physically-distinct phenomena. We emphasise that the former effect has been 
implicitly disregarded in previous KB studies that rely on particular quasi-particle ansaetze. 
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